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Abstract 

Bacteria and their viruses ("bacteriophages") coexist in natural environments forming 
Cln complex infection networks. Recent empirical findings suggest that phage-bacteria infec- 

•^T tion networks often possess a nested structure such that there is a hierarchical relation- 

ship among who can infect whom. Here we consider how nested infection networks may 
affect phage and bacteria dynamics using a multi-type Lotka-Volterra framework with 
cross-infection. Analysis of similar models have, in the past, assumed simpler interac- 
M-l tion structures as a first step towards tractability. We solve the proposed model, finding 

^H trade-off conditions on the life-history traits of both bacteria and viruses that allow co- 

existence in communities with nested infection networks. First, we find that bacterial 
growth rate should decrease with increasing defense against infection. Second, we find 



that the efficiency of viral infection should decrease with host range. Next, we estab- 



lish a relationship between relative densities and the curvature of life history trade-offs. 
We compare and contrast the current findings to the "Kill-the- Winner" model of multi- 
species phage-bacteria communities. Finally, we discuss a suite of testable hypotheses 
stemming from the current model concerning relationships between infection range, life 
history traits and coexistence in complex phage-bacteria communities. 



cr 



o 

f^i 1. Introduction 

m 

t-H Bacteria and their viral parasites, i.e., phages, are found in natural environments from 

oceans, soils to the human gut. There are an estimated 10 30 bacteria on Earth, with 

• i—i estimates of phages approximately 10-fold higher [TJI2]- Phages are not only abundant, 

but they are also key players in ecosystems. For example, phages can be responsible 
for a significant portion of microbial mortality, e.g., with estimates ranging from 20%- 
80% [SI HI [5] . These estimates of lysis are at the community scale. However, individual 
phages infect a subset of bacteria in a community. A growing number of empirical studies 
have begun to investigate the nature of cross-infections between phages and bacteria 
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[El |H1 H] ■ These studies have the potential to help identify the basis for phage- induced 
mortality, by delineating the specific phage types capable of infecting and lysing specific 
host types and, potentially, the taxonomic and biogeographic drivers of cross-infection 
|10j . Although predictive models of cross- infection remain elusive, it is evident that a 
single virus can infect multiple strains of a host [11] [12] , in some cases multiple species 
[6| [7] , and even hosts from different genera [13] . 

In an effort to identify cross-system trends and patterns in cross-infection, we (and 
collaborators) recently re-analyzed 38 phage-bacteria infection studies. This re-analysis 
identified a recurring pattern in these studies: viral strains have overlapping host ranges 
such that the overall network of infections is significantly nested [13] . In nested systems, 
there is a hierarchy for who can infect whom (see Figure 1A). In a perfectly nested system, 
the specialist virus can infect the most permissive host, the next most specialized virus 
infects the most permissive host and the second most permissive host, and so on (see 
Figure IB). Hence, the host that is most difficult to infect is infected only by the most 
generalist virus. In systems that are nested but not perfectly so, this nesting of infection 
ranges occurs more frequently than expected by chance. It is important to note that 
some of the nested phage-host infection networks re-analyzed in |14j are derived from 
experimental evolution studies in which the hosts and phages do not coexist at the 
same time point, but rather the nested relationship is only observed when performing 
cross-infection experiments between time points (e.g. [11). However, some networks 
are derived from ecological studies where samples are taken from the environment (e.g., 
[6l [Jj [12]). In such cases the finding of overlapping infection ranges poses a dilemma 
for understanding coexistence. In particular, how can a specialist virus coexist when 
it only infects a single host, and indeed the host that is most susceptible to infection? 
Further, how can a permissive host coexist with other more resistant hosts when there 
is a nested relationship to infectivity? The consequences of such interaction networks 
on the structure and dynamics of microbes and their viral parasites have not yet been 
established. Here, we integrate empirical observations of complex infection networks into 
ecological models of host-phage dynamics. 

Host-phage systems and their population dynamics have been studied mathematically 
for over 30 years [T51 fTB] . The earliest models were meant to facilitate understanding 
of relatively low-diversity chemostat experiments, often involving the cross-infection of a 
single phage type with a single host type. These studies often found that mutants could 
arise and so simple models were often extended to include two hosts and one phage, two 
phages and one host, and so on [T5], [TT]. However, host-phage models applied to natural 
environments require a greater diversity of bacteria and phage types. There are different 
types of approaches to integrate diversity into dynamic models. On the one hand, eco- 
evolutionary models have been developed that assume that new host and viral types 
can evolve as a result of ecological interactions. Despite their computational complexity, 
multiple examples of such eco-evolutionary approaches for phage-host dynamics are now 
available [TH1 HH HOI IUJ- Alternatively, ecological models in which a fixed diversity 
of types are included from the outset can be used to study how certain features in a 
community may help maintain diversity. Examples of such approaches include spatial 
multi-strain models [22] , models of competition between hosts possessing different types 
of immune systems [23] and models that incorporate higher trophic levels [53] . 

The most prominent theory of phage-bacteria ecological dynamics in multi-species 
communities is of the latter type and is known as the Kill-the- Winner model [2H [35] . In 
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this model, multiple species of bacteria and phage are considered, for which each virus 
is assumed to exclusively infect a single host type. The model's central conclusion is 
that each host type is controlled, in a top-down fashion, by a single viral type, with the 
exception of a single host type whose density is controlled by the total host biomass limit 
set by an additional generalist grazer. Moreover, the steady state densities of viral types 
are determined by relative differences in life history traits of the hosts (an issue we return 
to in the discussion). The model is meant to describe the dynamics of bacterial species, 
i.e., no strain level dynamics are considered, although its structure is general. 

Here we extend the basic framework of the Kill-the- Winner (KTW) model to in- 
corporate complex interaction networks and ask how nestedness mediates coexistence 
in multi-strain host-phage systems. We focus our attention on the idealized case of a 
perfectly nested interaction network, and show conditions on the life-history parameters 
necessary for coexistence. We find trade-off conditions necessary for coexistence and 
show that coexistence occurs even when the system is perturbed from the equilibrium. 
We examine the abundances of both hosts and viruses at steady state and their rela- 
tionship to life history traits and to infection range. Finally we examine the existence 
of coexistence equilibria in the general case where the infection matrix is not perfectly 
nested. We close by discussing the relevance of the current study to the KTW model, 
empirical efforts to link infection patterns with life history traits, and recent attempts to 
establish a link between network structure and biodiversity. 

2. Methods 

We model a system of n bacterial strains (hosts) and n viral strains (phages) . Bac- 
teria compete for implicit resources while viruses infect different subsets of the bacterial 
community. We denote the density of host i by Hi and the density of virus j by Vj . Our 
model of the ecological dynamics of the different host and viral strains is, 
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In the absence of viruses, hosts exhibit logistic growth with exponential growth rate r, 
and a community-wide carrying capacity K . To simplify the model we assume that for a 
given host, intrastrain and interstrain competition are the same. The parameters tfij and 
Pj are the adsorption rate and burst size (virion release per infection) of virus j (which 
we assume is independent of host i). Virus j decays outside the host at a rate m,j. M 
denotes the infection matrix, where My = 1 if virus j can infect host i and My = if 
virus j cannot infect host i. For analytical tractability, there are several characteristics 
of the life cycle of a viral infections that are not included in this model. For example, 
we assume there is no delay between infection and virion release. We also do not include 
the possibility of lysogeny, where viral genetic material is incorporated into the host 
chromosome and vertically transmitted to daughter cells for future activation and lysis. 

3 




21 19 17 15 13 11 9 7 5 3 1 



V irus 




Figure 1: (a) Infection network from an experimental study presenting a statistically nested pattern 
(original data from 1121 reanalyzed in |14p. The numbers identify different types of viruses and hosts, 
(b) Perfectly nested infection network. For a perfectly nested network, the numbers correspond to the 
rank (i.e number of interactions). White squares denote that a given virus can infect the host. 



Empirically measured phage-host infection networks are often statistically nested (for 
example, see Figure [la). Here, our main goal is to examine the ecological implications of 
nested interaction networks for host-phage dynamics. We start by studying the idealized 
case of a system where the infection network is perfectly nested; see Figure [TJa. In this 
representation of the network, viruses are numbered from the most specialist to the most 
generalist. We will refer to this number as the rank of a species. We see that virus 1 
infects only one strain of host whereas virus j infects j different host strains. Note that, 
in the nested case, virus j + 1 can infect all strains of hosts that virus j can infect plus an 
additional one. Similarly, hosts are ordered by the number of viral strains that can infect 
them. Host i can be infected by i different viral strains. Again, note that in the perfectly 
nested case, host i can be infected by all viruses that can infect host i + 1. We relax 
these simplifying assumptions in section |3.4| where we consider the case of interaction 
networks that are not perfectly nested (as is the case in Figure [1^,). 

Numerical simulations of the ecological dynamics were done using a Runge-Kutta 
method (ode45 in MATLAB 26 i). The parameters (life-history traits) used in all the 
simulations are shown in table |B.1| in |Appcndix B| These parameters were chosen from 
a baseline of biologically realistic values found in the literature [37J H5] . 



3. Results 

3.1. Equilibrium densities and conditions for coexistence 

We start by examining the equilibrium densities of the host and phage strains. When 
the 2n equations of system (nj) are simultaneously zero, we find the equilibrium densities 
corresponding to strain coexistence: 
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where the superscript * denotes equilibrium densities and hj 
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densities of the hosts are expressed exclusively in terms of the life-history traits of the 
viruses (hj). This is due to the top-down control of the hosts by the viruses. In Figure [2] 
we show graphically how to construct the host steady states. We see that virus 1, which 
exclusively infects host n, determines the equilibrium density of that host. The value of 
this density is the ratio hi = jy. This result is the same as in a system with just one 
host and one virus, and it can be interpreted as the host density necessary to support 
the virus infecting it (Figure pa). Virus 2, on the other hand, requires a host density 
of hi — J"? to persist in the system. This virus infects two hosts: one in common 
with virus 1 (host n), and an additional one (host n— 1). However, the density of host 
H n is already set by virus 1, so the density of host n — 1 is the difference between the 
required density for virus 2 (hi) and the density of virus 1 (hi) (Figure |2p) . Similarly, 
virus 3 infects one more type of host than virus 2, and thus, the density of the additional 
host (H n _2) is the difference between the density required by virus 3 to survive (/13) and 
the density set by virus 2 (hi); see Figure [2fc. The densities of the remaining hosts are 
determined in an analogous way, Eq pi). From this result we obtain that the total host 



biomass at steady state is J2 Hj = 



/3„0„ 



i.e. it is set by the most generalist virus. 
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Figure 2: Construction of the host steady states (H*). (a) The equilibrium density of host n, H*, is 
set by the aggregate life-history traits of virus 1, hi, which infects H n exclusively, (b) The equilibrium 
density of host n — 1, i?*_j, is set by the difference between /12 and hi. (c) The equilibrium density of 



Coexistence of all strains at positive densities requires three conditions to be met: 
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These conditions establish a connection between the life-history traits of the different 
viral and bacterial strains and the structure of the infection networks, specifically the 
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rank of the strains. Equation Q represents a tradeoff for the viruses and can also be 
written as h n > h n -\ > • • • > h% > h\. Larger values of hj imply a virus has a higher de- 
activation rate (higher rrij) and/or produces fewer viruses per infection and/or is worse 
at attaching to the host (smaller f3j(f>j)- Thus, the inequalities in equation Q provide 
the viral trade-offs necessary for coexistence: a virus with a broader host range has less 
advantageous life-history traits, as characterized by its ratio hj , compared to viruses with 
a narrower host range. Note that if two viral types infected a single host, the type with 
the lower ratio hj would out-compete the other. Here, coexistence is possible because the 
viral type with higher hj can infect more host strains than the viral types with lower hj . 
Equation (T5J) describes the trade-off between immunity and growth rate for the hosts. 
Specifically, coexistence is possible if a host that can be infected by more viral types 
has a higher growth rate compared to a host that can be infected by fewer viral types. 
Finally, equation (pi) specifies that the sum of the host densities at equilibrium needs to 
be less than the carrying capacity of the system. 

Figure [3] shows an example of the dynamics resulting from a system with 2 host 
and 2 viral strains which satisfies the conditions listed in equations Q, ^, and (J6J). 
In contrast, Figure [4] shows what happens when the conditions for coexistence are not 
satisfied. In this example, host 1 has a larger growth rate than host 2, which results in 
the extinction of host 2, and as a consequence, the extinction of virus 1. 

3.2. Community dynamics and invasion 

Conditions Q, |[5J), and ffib guarantee the existence of a coexistence equilibrium in 
system (HI . In this section we address two issues related to the coexistence between host 
and viral strains when those conditions are satisfied. First, via numerical simulations 
we investigate if the species densities tend to the coexistence equilibrium or if cyclic 
coexistence is also possible. Second, boundary equilibria where one or more of the host 
and viral strains are extinct also exist in system (nj) . We ask if those boundary equilibria 
are unstable with respect to invasion by the extinct host and viral strains when conditions 
Q, J5), and (jfj) arc satisfied. 

As seen in Figures [31 fsT and |6j when conditions Q, J5J), and ffity are satisfied, the 
host and viral densities can either tend to steady state or exhibit cyclic oscillations. In 
Figure [31 the coexistence equilibrium point is stable, and the values of the densities tend 
to the equilibrium values after transient oscillations. For a different set of life-history 
traits, the coexistence equilibrium point is no longer stable, but cyclic coexistence is still 
possible (Figure (5J) . Figure [6] is an example of cyclic coexistence for 5 species of bacteria 
and 5 species of virus that satisfy the trade-off conditions Q, (OH) , and ([6]) . In Figure |6l 
coexistence is still possible because the coexistence equilibrium exists, even though the 
time series may be irregular. Note that due to the dimension of the model we do not 
have a closed form solution for when cyclic dynamics arise in our system. 

In appendix A we show that conditions Q, (0)]), and (pi) imply that all boundary 
equilibrium points of system (nj) are unstable with respect to invasion by at least one 
host or viral strain that is absent from that subsystem. This implies that if the dynamics 
in that subsystem tend to the boundary equilibrium point, then that subsystem can be 
invaded by one or more of the extinct host or viral strains. A stronger conclusion can be 
reached if system (fTl) is permanent (i.e. densities are bounded above and, after some time, 
are bounded below by a finite value [29 ). Because the average long term invasability 
conditions along orbits in permanent Lotka-Volterra systems are equal to the invasability 
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Figure 3: Dynamics of 2 host types and 2 viral types when the infection network is perfectly nested and 
the life-history traits satisfy the conditions HI, tol, and JM for coexistence, (a) Growth rate, r;, as a 
function of host rank in accordance with condition J5J. (b) hi as a function of viral rank in accordance 
with condition HI. (c) Host densities as a function of time: both hosts coexist via a stable equilibrium, 
(d) Viral densities as a function of time: both viruses coexist via a stable equilibrium. 



conditions at equilibrium points [22], conditions Q, (p)]), and (p| imply that the extinct 
strains can always invade the subsystem when system (IT]) is permanent. It is an open 
question whether system (flj is permanent when conditions Q, (pi), and ffih are satisfied. 

3.3. Relationship between abundance and rank 

Another relevant question is the connection between rank and density. Is there a way 
to infer information about the infection network from measurements of density or vice- 
versa? We find that, with the exception of host n, host densities are determined by the 
difference between consecutive hi values (eq. (pi)). So, in our framework, measurements 
of density could inform us about the differences in the life-history traits of viruses with 
the most similar host range. If two viruses of consecutive rank are very similar (in terms 
of the aggregate life-history traits hi), then the corresponding host density is low. On 
the other hand, two very different values of hi for consecutive viruses imply high host 
density. The only exception to this analysis is the density of host n, which is determined 
uniquely by the ratio h\ of virus 1 and not by the difference between the ratios of viruses 
of consecutive rank. 

Figure [7] shows examples of the connection between the life-history traits of the viruses 
and the density of the hosts. In the special case where the trade-off curve for the hi has a 
curvature with a constant sign, the values of hi translate into a simple (monotonic) rule 
for the host densities as a function of rank. For example, when the trade-off between 
hi and viral rank is concave up, the H* increase with rank (Figure |7k). When the 
trade-off is concave down, the H* decrease with rank (Figure (7b) . In both cases, the 
density H* of host n can be an exception. Figure [Tj; shows the general case. The values 
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Figure 4: Dynamics of 2 host types and 2 viral types when the infection network is perfectly nested and 
the trade-off conditions are not met. One of the viral types does not survive because the life-history 
traits do not satisfy the conditions for coexistence, (a) Growth rate, r, as a function of host rank does 
not satisfy condition 151. (b) h as a function of rank, (c) Host density as a function of time: host 2 goes 
extinct, (d) Viral densities as a function of time: virus 1 goes extinct. 
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Figure 5: Cyclic coexistence of 2 host types and 2 viral types. The infection network is perfectly nested 
and the life-history traits satisfy the conditions for coexistence, (a) Growth rate, r^, as a function of host 
rank in accordance with condition (J5l. (b) hi as a function of viral rank in accordance with condition 
Q. (c) Host densities as a function of time: both hosts coexist via a stable equilibrium, (d) Viral 
densities as a function of time: both viruses are present via oscillations. 



3 

O 

rJ=i 0.5 



12 3 4 5 

Host rank 



1 b 




12 3 4 5 

Viral rank 




Time (hours) 



Figure 6: Cyclic coexistence of 5 types of hosts and 5 types of viruses. The infection network is perfectly 
nested and the life-history traits satisfy the conditions for coexistence, (a) Growth rate, r, as a function 
of host rank in accordance with condition 15b, (b) h as a function of viral rank in accordance with 
condition GJ, (c) host densities as a function of time: all of the hosts coexist, (d) viral densities as a 
function of time: all of the viruses coexist. 

of hi increase with rank following condition HI for coexistence, but the corresponding 
equilibrium densities need not be a monotonic function of rank. 

The viral densities depend, in part, on how different the host strains are. Specifi- 
cally, the equilibrium densities of the viruses depend on the differences between growth 
rates of host strains of consecutive rank. However, they don't depend uniquely on the 
traits of the host, they also depend on the adsorption rate of the focal virus (eq. (|3|). 
Thus, information about the viral densities need not translate directly into information 
concerning differences in the life-history traits of hosts. 



3-4- Deviation from a perfectly nested network and coexistence 

We now return to our assumption about interaction networks and consider systems 
that are not perfectly nested. The examination of imperfect nestedness is biologically 
relevant. For example, a re-analysis of 38 observational studies of phage-host interaction 
networks found numerous instances of elevated nestedness, all of which included some 
departures from perfect nestedness [14] . The ecological dynamics follow system (fil, 
however the different structures of the infection network will be reflected in the matrix 
M, which includes information of who can infect whom. 

The equilibrium densities can be expressed in a compact form as two matrix equations 
using the infection matrix M and its transpose M T : 
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Figure 7: Relationship between the viral tradeoff curve (left) and the corresponding host densities at 
steady state (right), (a) When the trade-off curve is concave up, host density decreases with rank, (b) 
When the trade-off curve is concave down, host density increases with rank, (c) In general, the host 
densities are set by the difference of consecutive hi and need not be a monotonic function of rank. 
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Figure 8: (a) Infection networks corresponding to an invertible infection matrix, (b) Time series from 
numerical simulation using the infection matrix in (a). 
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Here H* and V* are vectors of the equilibrium densities, ft, is a vector whose elements 
are the hi, and r is a vector whose elements are the r,. We also use the change of 
variable V*' = fiiV* . We are interested in the solutions of system Q that are positive. 
We consider two cases: the infection matrix is invertible and the infection matrix is 
singular. An invertible infection matrix can be interpreted biologically as each viral 
strain having an unique niche, where the niche is defined by the host range. A singular 
infection network can be interpreted as the existence of niche overlap between two or 
more viral strains. 

First consider the case where M is invertible, i.e., there is niche differentiation. In this 
case, there exists a unique solution to system n\, and therefore a unique set of coexistence 
equilibrium densities. The equilibrium densities are expressed in terms of the life-history 
traits of the different viruses and bacteria and will, in general, involve differences between 
traits or combinations of traits (e.g. hi). Therefore, for every invertible infection network 
there exists a series of inequalities that are necessary to guarantee positive equilibrium 
densities (analogous to conditions Q, (pjl), and Q for the perfectly nested case). In 
Figure [8k we show an example of an infection network that corresponds to an invertible 
infection matrix that is not perfectly nested. We also show time series for a system 
with that specific interaction matrix where all the strains are present in the community 
(Figure pb). The relationship between network structure and equilibrium densities may 
be suitable for further numerical analysis. 

If the matrix M is not invertible, i.e. there is niche overlap between viral strains, 
then there are two ways in which coexistence can occur. The first case occurs when the 
life-history traits of some viral strains are effectively equal. For example, in a system 
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with two viral strains and one host strain, the viral strains can coexist if they have the 
same hi values. Biologically, this case is unlikely to occur given that we are describing 
strains in terms of function, and in this case, coexistence is possible only when strains 
are functionally identical. 

The second case in which coexistence is possible occurs if we relax assumptions about 
the adsorption rate and the burst size for viral strains. In our formulation, a specific viral 
type infects all host strains with the same adsorption rate and burst size. However, if 
viruses exploit host strains at different rates, i.e., adsorption rate and burst size depend 
on both viral and host type, then coexistence is possible even with complete niche overlap. 
In this case, the steady states are solutions of two analogous systems of equations where 
the infection matrix M is a weighted matrix whose entries correspond to the adsorption 
rate and burst size for each host-viral interaction. A limited number of studies have 
shown that viral infection rates can differ significantly between host strains (e.g., see 
[13]). Hence, biologically, the use of quantitative information for host-phage infection 
assays represents an important target for future analysis and may shed light on the 
drivers of coexistence in natural communities. 

4. Discussion 

We studied the ecological dynamics of phages and hosts using a Lotka-Volterra frame- 
work that incorporated complex cross-infection networks. We found that coexistence is 
possible even when viruses exploit overlapping ranges of host. In the case of a perfectly 
nested infection network, Figure [lb, we found trade-offs for both the hosts and the viral 
strains that are necessary to allow the coexistence of all the species. The trade-off for 
hosts implies that host growth rate decreases with defense (the larger the growth rate, 
the larger the viral range). The trade-off for viruses implies that viruses must be less 
efficient at utilizing host resources as they increase their host range. We also showed 
that the densities of the host strains at equilibrium are determined by how different the 
viruses are in terms of their aggregate life-history traits. 

The idealized case of a perfectly nested network is mathematically tractable and may 
help to identify potential principles underlying coexistence in host-phage systems. On 
the other hand, real infection networks are rarely perfectly nested (see Figure 1 for an 
example, and the re-analysis of Flores et al [2]). We showed that coexistence is possible 
in the current framework if there exists a partitioning of niches. In those cases, the 
resulting tradeoffs between infection range and life history traits are not easily presented 
in a general way because each infection network results in different trade-off conditions. 
This, in turn, makes it difficult to study the stability for the general case. Nevertheless, we 
showed numerical examples of coexistence away from equilibrium for a infection matrix 
that is not perfectly nested. 

The current model is similar in spirit to the Kill-the- Winner (KTW) model [21] • The 
KTW model also proposed mechanisms by which bacteria could coexist, stabilized by 
the presence of viruses. However, in the KTW model, infections are one-to-one, meaning 
that each bacterial strain can only be infected by one virus and likewise, each virus can 
only infect one bacteria. This would be equivalent to a diagonal infection network in our 
representation. In the KTW model, coexistence of different types of bacteria is achieved 
through what is known as a "killing the winner" mechanism, where "coexistence among 
bacteria is ensured by host-specific viruses that prevent the best bacterial competitors 
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from building up" |24) . Hence, viruses enable coexistence in a system that would, in 
their absence, lead to a diversity collapse. The principle of top-down control also applies 
to the model presented here, despite the fact that we considered complex interaction 
networks with the possibility of overlapping host range. 

However, the current model makes predictions not found in the KTW model. First, 
KTW assumes each bacterial strain is infected by a single viral type and so an ordering 
of bacterial growth rates is not necessary for coexistence. In the current model, an or- 
dering of bacterial growth rates with respect to susceptibility to infection is required for 
coexistence and represents a prediction of the model. Second, whereas KTW assumes 
no difference in viral life history traits, we again predict that differences in life history 
traits are required for coexistence when viruses differ in their host range. Hence, alto- 
gether we predict that there should be an entanglement between network structure and 
bacterial/viral life history traits. 

In light of the differences in predictions between our models and the KTW model, we 
suggest that studies of cross-infection from the enviornment should move beyond quali- 
tative analysis (i.e., whether or not a phage infects a bacteria) to quantitative analysis 
(i.e., the lysis rate of bacteria by phages). Doing so would help identify costs of resistance 
and infectivity in natural populations as well as help understand the relative importance 
of network structure and trade-offs in shaping the structure of interactions in natural 
communities. We note that extensions of the current model will be required to consider 
the type of trade-offs that may coincide with findings of imperfectly nestedness infection 
networks. 

In the past, the relationship between infection range and physiological costs has been 
studied experimentally, often in model organisms. For example, Lenski studied 20 E. 
coli mutant strains resistant to phage T4 [30 . All 20 strains were significantly less fit 
that the parental strain as measured via direct competition. Similarly, Bohannan and 
colleagues found strong growth costs to resistance for E. coli strains resistant to phage 
A and to phage T4 [31j . The hypothesized mechanism is that resistance often requires 
modification to surface receptors that may alter the rate and effectiveness of nutrient 
uptake [3"2"] . 

However, resistance mechanisms need not depend on surface changes and, moreover, 
in practice such trade-offs are not found universally. Indeed, two recent studies of putative 
trade-offs found more equivocal evidence. In one instance, Lennon et al |33j studied 
Synechococcus hosts and associated myoviruses. They considered 22 different bacterial 
strains selected for viral resistance, based on four different ancestral strains. 11 out of 
the 22 strains showed a fitness cost compared to its ancestral strain either in the form of 
its maximum growth rate from growth curve data or through direct competition assays. 
Hence, 11 of the 22 strains showed no fitness cost to resistance. Similarly, Avrani et al [34] 
studied the relationship between Prochlorococcus hosts and associated podoviruses. They 
found that 11 of 23 mutant strains that had evolved phage resistance grew significantly 
more slowly than did the ancestral wildtype. Hence, 12 of the 23 mutants strains did not 
show a growth cost to resistance. Instead, they found an alternative type of trade-off in 
which resistant mutants were susceptible to infection and lysis at a higher rate by other 
podoviruses and myoviruses with which they had not coevolved [34 . Hence, it may be 
that the costs of resistance can not be fully understood without taking the community 
into account, as was considered in the model presented here. 

There have also been a limited number of studies linking host range expansion to 
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fitness cost within phages. For example, Duffy and colleagues showed that the fitness of 
evolved RNA-based phages, when grown on the ancestral host, declined significantly in 
28 of 30 phages that had evolved an expanded host range [55] ■ However, 2 of the phages 
had no significant fitness cost [35] . Given the diversity of hosts and viruses in natural 
systems, it seems important to extend prior assays relating one-step range expansions 
to fitness. In particular, our work and the above suggest it is important to measure 
both host range and physiological costs in relationship to a broad range of potentially 
interacting strains. 

The integration of life history trait variation into studies of community composition 
requires taking into account the coevolution of phages and bacteria. Previous studies have 
explored the coevolution of phage-host theoretically [HI |36l HH [32 [201 E] ■ Experimental 
studies of co-evolving phage-host systems suggest that infection networks are not static, 
but rather are dynamic and reflect the changing identities of strains in a population 
[581 1551 HP] . Phage-bacteria coevolution is often described in terms of an arms race where 
the species' traits escalate over time [53 [52]. In such models it is common for better 
defended host types and phage types with better offenses to replace the respective resident 
type. However, coexistence of types is also possible if disruptive selection or evolutionary 
branching occurs. We have shown that ecological coexistence between multiple phage 
and bacterial types is possible if a series of trade-offs are satisfied (see Equations Q- 
([6])). Thus, if there is an appropriate trade-off between host defense and growth rate (e.g., 
resistance to infection) and an appropriate trade-off between phage replication and host 
range (e.g., ability to bypass host defense), evolutionary branching could yield phage- 
bacteria communities with nested patterns of infection. Thus, an arms race between 
phage and bacteria could result in the coexistence of one phage and bacteria pair that 
evolves over time or the coexistence of multiple phage and bacteria types that lie along 
their respective trade-offs. 

We also point out that coevolution also occurs over spatially extended domains. Spa- 
tial structure is thought to stabilize diverse interactions amongst phage and bacteria 
(e.g. [JTJ (25J [55]). Moreover, geographic structure can play a key role in affecting the 
outcome of coevolution [42]. Further work is warranted to quantify how structure in 
infection networks are driven by and act as drivers of the spatial distributions of diverse 
communities of phages and bacteria [43j HU HS] . 

Altogether, the current study predicts that phage-host coexistence in a given system 
depends on both network structure (e.g., nestedness) and life history traits. Other types 
of infection networks are possible including modular |14] and multi-scale [431 , suggesting 
the need for further investigations into the relationship between cross-infection and life 
history traits. Extending the current model framework to include evolution, spatial dy- 
namics, and quantitative information on infection will yield increased opportunities for 
theories of phage-host interactions to deepen our understanding of the mechanisms un- 
derlying the formation of complex phage-host interactions networks and their ecosystem- 
level consequences. 
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Appendix A. Stability of the boundary fixed points for phage-host dynamics 
with a perfectly nested infection matrix 

In this section we show that if conditions Q, (pH), and (pf from the main text are 
satisfied, then all boundary equilibria of the system are unstable. As discussed in the 
main text, this implies that if the dynamics in the boundary subsystem tend to the 
boundary equilibrium, then the subsystem can be invaded by at least one of the extinct 
host or viral strains. If trajectories in the subsystem do not converge to the boundary 
equilibrium point, then our conclusions about invasability only apply if system (fTl) is 
permanent (i.e. densities are bounded above and, after some time, are bounded below 
by a finite value, |29jV When system (fTl) is permanent, invasion at the equilibrium 
implies invasion along any orbit (over infinite time). This is because the average long 
term invasability conditions along orbits in permanent Lotka-Volterra systems are equal 
to the invasability conditions at equilibrium points |29| . We are not aware of any proof 
that our system is or is not permanent. 

For notational convenience, we make two changes to the way the system is presented. 
First, we don't use the matrix Mjj to establish who can infect who. Instead we change 
the range of sums in accord with the nested pattern. Second, we reverse the order of the 
hosts (Hi — ¥ H n _ i+ i). In this numbering, host 1 can be infected by all viral strains, and 
host n can be infected by a single viral strain (Figure A.l). In this new notation the 
system takes the form: 



I E H A 




-HiJ^tM, 



(A.l) 



(A.2) 



In this notation, conditions Q, ^, and ^ from the main text become: 



m 


+1 


> 


m-i 


<i>i+i 


A+i 


M 


n+i 


< n 






m n 


^ 


r s 





(A.3) 
(A.4) 
(A.5) 



Note that, as a consequence of reversing the order of the hosts, condition (A.4) on the 
growth rates (r,) is the reverse of condition ^. 

We will show that all possible boundary fixed points are unstable with respect to 
invasion. Let X = (Hi, H 2 , ■ . ■ , H n ,V\, V2, ■ ■ ■ , V n ) denote the vector of all host and viral 
densities and let X* = (H % , J?| > • • • > H*, V* , V 2 * > • • • » V£) denote a fixed point of system 



( A.l I and (A. 2), where H* is the equilibrium density of host i and V* is the equilibrium 
density of virus i. We will denote the Jacobian by J and the eigenvalues by A. We will 
consider 4 properties of fixed points, such that every possible boundary fixed point has 
at least one (and possibly more) of these properties. We will show that any fixed points 
with any of these properties is unstable with respect to invasion by at least one of the 
host or viral strains. The 4 properties of fixed points are: 

1. X* where H? = V i E [l,n]. 

2. X* where H{ = H* 2 = H 3 * = • ■ ■ = H * k = with k < n and H* k+1 ^ 0. 

3. X* where H* = H* +1 = • • • = H* = for i^ 1 and H t ^ ^ 0. 

4. X* where H* = H* +1 = ■ ■ ■ = H* k = for i > 1, k < n, and H k+1 ^ 0. 

The first property corresponds to fixed points where none of the hosts are present. 
The second property corresponds to fixed points where hosts with low defense level are 
not present (i.e. hosts that can be infected by all or many viral types are not present). 
The third one corresponds to fixed points where hosts of high defense are not present 
(i.e. hosts that can be infected by one or only a few viral types are not present). The 
last property corresponds to fixed points where host of intermediate defense level are not 
present. 

We will prove the instability of fixed points with properties 1,2,3, and 4 in lemmas 3, 
4,5 and 6 respectively. Before that we will prove two more lemmas. 

Lemma 1. Let X be a fixed point, then: 

a) If H* = 0, then ij^ „» is an eigenvalue of the Jacobian evaluated at X* . 

b) If V* = 0, then ^y:\ x , is an eigenvalue of the Jacobian evaluated at X* . 

Proof, a) 

Assume H* = 0. We can write the equation describing the dynamics of Hi as: 

it W 

^ = J?i/(X), (A.6) 

where /(X*) is a function of all the densities. Then, for any density X 7^ Hi, 

= 0. (A.7) 



m 

dX 



,-W) 



1G 



x* 



Hence, the Jacobian will have the form: 
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/ am 
' am 








am 

dH„ 


dHi 


... 





dm 
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dV„ 

\ dm 








dV n 
8H n 


dV ± 



em \ 

dV n \ 



(A.5 



dV n , 
dV n / 



X* 



Note that all the elements in the ith row are zero, except for the diagonal term, which 

jristl 

dHi 



has the form §g*| x „. Thus, the characteristic equation can be written as 



det( J - IX) = 



dH,, 



where P(X) is a polynomial A. Hence, W^ 

b) The proof for the case where V* 
H* = 0. 



- A P(A) = 0, (A.9) 



is an eigenvalue of the Jacobian. 

follows the same steps as proof a) for 

□ 



Lemma 2. Let X* be a fixed point where V* = for all i and H* ^ for at least one 
i, then X* is unstable with respect to invasion by at least one virus. 



Proof. Assume H* ^= 0. The equilibrium condition for Hi is 



Ti\ 1 



K 




0. 



(A.10) 



So, EH*=K. 

Now, V* = 0, so using lemma [T] we know that one eigenvalue of the Jacobian is: 



A = 



dv: 



0V n 



x* 



= 0n/3n ^2 H* - m n = <j> n p n K - m n > 0, 



(A.11) 



where the inequality results from the condition of coexistence (A.5). Equation (A. 11) 



implies that the fixed point X is unstable with respect to invasion by virus n. □ 

Now, we will show that fixed points with any of these 4 properties are unstable. 

Lemma 3. Let X* be a fixed point where H* — for all i e [l,7l], then X* is unstable 
with respect to invasion of any of the hosts. 



Proof. Assume H* = for all types of hosts in the system, then from equation ( A.l ) it 
follows that also V* = for all types of viruses at equilibrium. 
Using lemma [T] 



A 



dHj 

dH,, 



n>0 



X' 



is an eigenvalue. Thus X* is unstable with respect to invasion of any of the hosts 
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(A.12) 
□ 



Lemma 4. Let X* be a fixed point where H{ = H^ = ff| = • • • = H%. = with k < n 
and H^ +1 ^ 0, then X* is unstable with respect to invasion of hosts 1 through k . 

Proof. Assume H± = H^ = H£ = ■ ■ ■ = H£ = with k < n. The equilibrium condition 
for Vi for 1 < i < k is: 



,ft Y, H J -mi]Vi = 0. 



(A.13) 



Note that V\, V~2, ..., Vk only infect hosts 1 through k (the upper limit of the sum is i and 
1 < i < k ). Thus if H* — for the first k viruses, then it must be the case that V* = 
for the first k viruses, i.e 



VI =V£ = V£ = ■■■ = V k * = 0. 
From lemma [T] the eigenvalues for the first k hosts have the form 



X = 



dH z 



A'* 



>-W)-± 



W 



(A.14) 



(A.15) 



j=k+i 



Note the lower limit in the sum as a result of equation (A.14). 



Now, assume Hk+x 7^ 0. From the equilibrium condition for H^+i we get: 



j=k+i 



r k + 



^-W 



(A.16) 



If V* = for all i in the sum, then V* = for all the viruses in the system. This case is 
already covered in lemma [2] 

On the other hand if V* 7^ for at least one virus V* we get, 



'-¥ 



>0. 



Plugging the result (A.16) in equation (A.15), we get: 

Eh; 



A = {n - r k+1 ) 1 



K 



>0, 



(A.17) 



(A.18) 



where the inequality comes from the inequality (A.17) and the condition for coexistence 



(|A.4|). Equation (|A.18|) implies that the fixed point X* is unstable with respect to 

□ 



invasion by host i where 1 < i < k. 



Lemma 5. Let X* be a fixed point where H* = H* +1 = ■ ■ ■ = H* = for i =/= 1 and 
H*_ 1 =/= 0, then X* is unstable with respect to invasion of at least one species of host or 
virus. 
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Proof. Assume H* = H* +1 
through n are: 



H* — 0. The equilibrium condition for viruses i — 1 



E^ 



mj_i 



8-1 



v.' =1 



0zft 



K-i = o 



K* =0 



(A.19) 



E^; 



^nA, 



K? = o. 



Note that all the sums in (A.19 1 are the same. Because -Jpi- ^ npg- for all £ ^ p in 



the system, the equations in (A.19) can only be satisfied if V£ ^ for at most one 
k e [i—l,ri\. 

To see that X* is unstable, consider three possible cases: 

Case 1: X* where V fe * = V k G [i - l,n]. 

Case 2: X* where V* ^ and V fc * = V fe € [i - 1, n - 1]. 

Case 3: X* where V fc * 7^ for one fee [i — 1, n — 1] and Vg = for € G [i — 1, n], k ^ t. 

Cased) 

Assume V fe * = for all fe G [* — 1, n]. In this case V fe * = for all k G [l,n]. Indeed, 
the equilibrium condition for Hi-\ implies: 



4-W 




Thus, 



E h ;= k - 



(A.20) 



(A.21) 



From equation (A.l) we see that equation (A.21 1 implies V* = for all i. This is the 
case covered in lemma (J2J). 

Case[U 

Assume V* 7^ and V k * = for all fc G [i — l,n — 1]. From lemma 1 we know that, 
because Vk = for i — 1 < k < n the corresponding eigenvalues have the form: 



x _on 



dV k 



i-\ 



WkJ2 H l 



TOfe. 



X * 3=1 

From the equilibrium condition for V n , we get: 



E^; 
3=1 



(A.22) 



(A.23) 
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Note that the limits of the sums in (A.22) and (A.23) are the same as a result of the H 



that are zero, substituting the last expression into (A.22 1 

m 



A = 0fc/3fc 



4>nfin 



m k > 0, 



(A.24) 



where the inequality results from (A. 3 1, and consequently the steady state X* is unstable 



with respect to invasion of any of the viruses i — 1 through n — 1. 

Case [3] 

V fe * 7^ for one fee [i-l,n-l] and Ve = for I e [i - 1, n], ft 7^ 1 The only virus 
capable of infecting host n is virus n. Using lemma 1 we get that one eigenvalue is given 
by: 



A = 



8H n 



dH n 



x* 



1- 



Efi; 



K 



(A.25) 



where we used V* — 0. At this fixed point the only virus infecting host i — 1 is virus ft, 
so the equilibrium condition becomes: 



n-i 1 



K 



Thus, 



and 



1- 



Eh; 



K 



4> k v k * = 0. 



>o 



K 



>0. 



So, the fixed point X* is unstable with respect to invasion by host n. 



(A.26) 

(A.27) 

(A.28) 
□ 



Lemma 6. Let X* be a fixed point where H* = H* +l = ■ ■ ■ = H^ = for i =/= 1, k < n, 
and Hk+i 7^ 0, then X* is unstable with respect to invasion of at least one species of host 
or virus. 

Proof: 

Similarly to lemma [5] the equilibrium condition for viruses i — \ through ft are: 



20 



i-1 

rrii-i 



E^-Srh-^ ( A - 29 ) 



j =i 



4>i-iPi-i 



£*-W*- 



To simultaneously satisfy all of the equations at most one Vf ^ for £ £ [i — 1, k] ■ 
To show that X* is unstable, it is convenient to study two different cases. 

Case 1: X* where F fc * = 
Case 2: X* where V fe * ^ . 

Cased) 

Assume Vfe = 0. Using lemma 1 we know that an eigenvalue of the Jacobian is: 

fe \ / j=k+l 

The equilibrium condition for Hk+i yields: 

rk+1 V'K )- £ W = °- ( A - 31 ) 

v 7 j=fe+i 

Using this result the eigenvalue becomes: 



,)d-5J£ 



A = (r fc - r fe+1 ) ( 1 - ^- ) . (A.32) 



And as a consequence of condition (A.3), the eigenvalue is positive and the fixed point 
is unstable with respect to invasion by host k. 

CaseE] 

Assume Vfc 7^ 0. The equilibrium condition for Vfc is: 

i-l 

Af3 k J2 H j" m k = 0. (A.33) 

3=1 



So, 
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From Lemma 1 follows that for any i < £ < k, there exists an eigenvalue of the form: 

i-l 



A 



av; 



dV e 



x> 



&Pij2 H i 



me. 



Substituting equation (A.33), we get: 

8V, 



A = 



dV f 



j'=i 



m k 



= (f>tfa , , 

X- 9kPk 



mi > 0, 



(A.35) 



(A.36) 



where the inequality is a consequence of condition (A. 3 1. Thus, the fixed point is unstable 



with respect to invasion by viruses i — 1 through k — 1 
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Appendix B. Parameters values used in the numerical simulations 







r (h-^) 


K (colls/ml) 


mth- 1 ) 


(ml/ (cells xh)) 


(3 (No. of viruses) 


Fig(3] 


Hi 


0.1760 


10'' 










H 2 


0.9914 


10' 










Vi 






0.0681 


5 x 10~ y 


24 




V 2 






0.1047 


5 x 10~ y 


10 


Fig 


4 




Hi 


1.0500 


10 s 










H 2 


0.300 


10 s 










Vi 






0.0100 


5 x 10~ a 


100 




V 2 






0.0650 


5 x 10 -9 


100 


Figg 


Hi 


0.3000 


10° 










H 2 


1.0500 


10 s 










Vi 






0.0100 


5 x 10~ u 


100 




v 2 






0.0650 


5 x 10~ u 


100 


Fig 


6 




Hi 


0.3162 


10' 










H 2 


0.6325 


10' 










H 3 


0.8367 


10' 










Hi 


1.1000 


10' 










H 5 


1.1400 


10' 










Vi 






0.0090 


5 x 10 _a 


20 




v 2 






0.0270 


5 x 10 _a 


40 




v 3 






0.0548 


5 x 10 _a 


60 




Vi 






0.0922 


5 x 10 _a 


80 




v 5 






0.1392 


5 x 10 _a 


100 


Fig|J> 


Hi 


0.4472 


10' 










H 2 


0.7583 


10' 










H 3 


0.9747 


10' 










Hi 


1.1511 


10' 










Vi 






0.0100 


5 x 10~ u 


20 




V 2 






0.0375 


5 x 10~ u 


46 




v 3 






0.0650 


5 x 10~ u 


70 




Vi 






0.0925 


5 x 10~ u 


96 



Table B.l: Values of the life-history traits used in the numerical simulations. Values based on |27l I28j 
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